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A Note on a Three Variables Analogue of Bessel Polynomials 

Bhagwat Swaroop Sharma 

Abstract: The present paper deals with a study of a three variables analogue of Bessel polynomials. 
Certain representations, a Schlafli 's contour integral, a fractional integral, Laplace transformations, some 
generating functions and double and triple generating functions have been obtained. 



I. Introduction 

In 1949 Krall and Frink [12] initiated a study of simple Bessel polynomial 

X 



Y n (x) = 2 F D 



n,l + n;-; 



and generalized Bessel polynomial 



Y n (a,b,x)= 2 F D 



n, a - 1 + n: 



x 
b 



(1.1) 



(1.2) 



These polynomials were introduced by them in connection with the solution of the wave equation in 
spherical coordinates. They are the polynomial solutions of the differential equation. 

x 2 y" (x) + (ax + b) y' (x) = n (n + a - 1) y (x) (1.3) 
where n is a positive integer and a and b are arbitrary parameters. These polynomials are orthogonal on the unit 
circle with respect to the weight function 

f 2 ^ 



1 



p(x,a) = — Yj 



r(a) 



27ii n=0 T(a + n - 1) 



(1.4) 



Several authors including Agarwal [1], Al-Salam [2], Brafman [3], Burchnall [4], Carlitz [5], 
Chatterjea [6], Dickinson [7], Eweida [9], Grosswald [10], Rainville [15] and Toscano [19] have contributed to 
the study of the Bessel polynomials. 

Recently in the year 2000, Khan and Ahmad [11] studied two variables analogue Y^"'' 3 ' (x, y) of the 
Bessel polynomials Y^ (x) defined by 



Y< a) « = 2 F 0 



n, a + n + 1; 



(1-5) 



The aim of the present paper is to introduce a three variables analogue Y^"'^ (x, y ,z;a, b, c) of (1.2) and to 
obtain certain results involving the three variables Bessel polynomial Y i | a '^' y ' ) (x, y ,z;a, b, c) . 



(a,p,y) 



II. The Polynomials 

Y n (a ' P ' Y) (x, y,z;a,b,c): The Bessel polynomial of three variables Y^'"" (x, y,z;a, b,c) is defined as 

follows: 
Y n (oJW (x,y,z;a,b,c) 

=Z £ n £ ( -~ n ^ +s+ ^ a+n+1 ^ p+n+1 ^ Y+n+1 - >i - —— -- (2D 

r =o s =o j=o r ! s ! j ! 

For z = 0, a = b = 2, (2.1) reduces to the two variables analogue Y^ a ' P) (x, y) of Bessel polynomials (1.5) as 
given below : 

Yf^ (x, y,0;2,2,c) = (x, y) (2.2) 

Similarly 

Y ( n aM (x, 0, z; 2, b,2) = Y n (a ^ (x, z) (2.3) 





x^ 
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Yf p " ) (0,y,z;a,2,2) = Y^(y,z) 
Also for a = - n - 1 , and b = c = 2 

Y(-n-i>P>y) (x , y , Z ; a ,2,2) = Yf^(y,z) 

Similarly, 

Y f- n -^(x,y,z;2,b,2) = (x, z) 

Y (a,p,-n-i) (X) y;Z;2 , 2 ,c) - Y n (a ' P) (x, y) 

where 

Y (a,p )(x?y) ^^ g (-n) r+s (a + n + l) s ((3 + n + l) r 
r =o s=o r ! s ! 

Also, for y = z = 0, a = 2, (2.1) reduces to the Bessel polynomials Y^ a) (x) as given below: 

¥<***> (s, 0,0; 2, b,c) = Y,: a) (x) 

where Y^ a) (x) is defined by (1.5). 
Similarly 

Yf M (0,y,0;a,2,c) = Yf (y) 

^^(aazA^c) = y„ w (z) 

Also, for p = y = - n - 1, a = 2, we have 

Yf-- n - 1 '- n - 1) (x,y,z;2,b,c) = Y { n a) (x) 

Similarly 

Y (- n - 1 -P-- n - 1 >(x,y,z;a,2,c) = Y<» (y) 

Y^ n - 1 -- n - 1 ' 1 ' ) (x,y^;a,b,2) - Y«(z) 



(2.4) 
(2.5) 

(2.6) 
(2.7) 

(2.8) 
(2.9) 



(2.10) 
(2.11) 

(2.12) 

(2.13) 
(2.14) 



y+n 



III. Integral Representations 

It is easy to show that the polynomial Y^ a '^' Y) (x, y,z;a,b,c) has the following integral representations: 
1 c 

\ r r r u <*+" v p +n w 

r(a + n + l)r(|3 + n + l)r(Y + n + l) Jo Jo Jo 

= Y (»JW (x,y,z;a,b,c) 
For z = 0, a = b = 2, (3. 1) reduces to 

1 < 



, xu yv zw 
1 + — + — + — ] e 
a b c 



uvw dudvdw 
(3.1) 



JtCO poo 
u 
o Jo 



a+n yP+ n 



, xu yv 

1 + — + — 
2 2 



V 1 



r(a + n + l)r((3 + n + l) 

= Y n (a ' P) (x y) 

a result due to Khan and Ahmad [11]. 

For y = z = 0, a replaced by a - 2 and a replaced by b, (3.1) becomes 

e ' dt 



1 + * 

V by 



(3.2) 



(3.3) 



a result due to Agarwal [1]. 

£ \l £ x a (r-xr 1 yP(s-yr 1 z^(t-zr 1 Y^)(x,y,z;a,b,c)dxdydz 

a+n p+n t Y+n /t-/ Y)3 

_ r s t {L{n)j ( a - n ,p- n , Y - n ) / , \ 

-7 7\ 7z — 7\~ x n ^r,s,t,a, D,cj (3.4) 

(a + l) n (p + l) n (y + l) n 

Jo Jo Jo ^l 1 " 11 )" v l4 (l-v) M w^'(l-wf 1 Yi a * T) (xu, yvzw;a,b,c)dudvdw 
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_ r(5)r(^)r(T 1 )r(^)rfe)r(v) F(3) 
r(s+^)r(Ti+n)rO;+v) 



-n::-;-;-:a+n+l,8;P + n+l,Ti;y+n+l,^; x y z 

-::-;-;-: 8 + A, ; r\ + [i ; £, + v ; a' b' c 

(3.5) 

where F <3) [ ] is in the form of a general triple hypergeometric series F <3) [x, y, z] 
(cf. Srivastava [18],p. 428). 

f £ Jo u 5 ^ 1 " 11 )^ 1 v^'(l-v) M w^(l-w) v l Y i ! aPy) (x(l-u),yvzw;a,b,c)dudvdw 



_ r(5)r(x)r(T 1 )ryrfe)r(v) F(3) 
r(5+^)r(T 1 +^)rfe+v) 



i::-;-;-:a+n + l,X;P+n + l,Ti;y+n + l,^; x y z 

::-;-;-: 8 + X ; r\ + n ; £, + v ; a' b' c 

(3.6) 

Jo Jo Jo ""(i-ur 1 v^l-v)^ 1 w^l-w)^ 1 Y n (aPy) (xu,y(l-v),zw;a,b,c)dudvdw 



r(8+x)r(T,+n)rfe+v) 



7(3) 



-n::-;-;-:a + n + l,5;p + n + l,|u;Y + n + l,^; x y z 

-::-;-;-: 8 + A, ; r| + ju ; £, + v ; a' b' c 

(3.7) 



Jo Jo Jo u 5 ^ 1 " 11 )^ 1 v^'(l-v) M w^(l-w) v l Y 1 ! aPy) (xu,yv,z(l-w);a,b,c)dudvdw 

_ r(5)r(x)r(T 1 )r(^)rfe)r(v) F(: 



■n::-;-;-:a+n+l,8;P+n+l,Ti;y + n+l, v; x y z 

-::-;-;-: 8 + X ; r\ + n ; ^ + v ; a' b' c 

(3.8) 



f f f u^fl-uf" 1 v^l-v)^ 1 w^(l-w) v "' Y n (a ' p ' y) (x(l-u),y(l-v),zw;a,b,c)dudvdw 

JO JO JO 



_ r(8)r(x)r(T 1 )ryrfe)r(v) F(3) 
r(8+^)r(Ti+^)r(^+v) 



i::-;-;-:a + n + l, A,;P + n + l,|u;Y + n + l,^; x y z 

::-;-;-: 8 + X ; r| + ju. ; <; + v ; a' b' c 

(3.9) 

Jo Jo Jo u6_1 ( 1 - u ) X_1 v' H (l-v) M w^fl-w)^ 1 Y,! aPy) (x(l-u),yv,z(l-w);a,b,c)dudvdw 
_ r(8)r(x)r(ii)r(^)rfe)r(v) p(3) r-n::-;-;-:a+n+lA;P+n+l,Ti;y+n+l,v;_x _y _z 
r(8 + X)r(r\ + |u)r(^ + v) |_ -::-;-;-: 8 + X ; r) + |u ; t, + v ; a' b' c 

(3.10) 

f f f u 51 (l-u)" 1 v n l (l-v) M w^(l-w) v l Y I J aPy) (xu,y(l-v),z(l-w);a,b,c)dudvdw 

JO JO JO 

i::-;-;-:a+n+l,8;P+n+l,|4,;y + n + l,v; x y z 
::-;-;-: 8 + A, ; r\ + \x ; ^ + v ; a' b' c 

(3.11) 



r(8M^MnMQr(v) 
r(8+^)r(T 1 +^)rfe+v) 



p(3) 



_[ J J u 8 " 1 (l - uf~' V th (l - v) M w 5-1 (l - w) v -' Y n (a P y) (x(l - u), y(l - v), z(l - w); a, b, c) du dv dw 

_ r(8)r(x)r(T 1 )r(^)rfe)r(v) F(3) 
r(8+x)r(T 1 +^)rfe+v) 



-n :: 



;-:a + n + l,A,;P + n + l,|u;y + n + l,v; x y z 

-: 8 + A, ; r\ + n ; <; + v ; a' b' c 

(3.12) 

£ £ £ u s -Xl-u) x -\^Xl-v) M w^ 1 (l-wr 1 Y^(xvw,yuw,zuv;a,b,c)dudvdw 



_ r(8)r(^)r(T 1 )r( H )rfe)r(v) F(3) 
r(8+^+^+v) 



-n::8 ; r) ; ^:a + n + l;P + n + l;y + n + l; x y z 

-:: 8 + X ; rj + fa. ;^+v : -; -; - ; a' b' c 

(3.13) 
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( ( ( u 51 (l-u)" " 1 v ,1 - | (l-v) M w^ 1 (l-w) v " 1 Y< a ' Py) (xv(l-w),yw(l-u),zu(l-v);a,b,c)dudvdw 

JO JO JO 



_ r(5)r(x)r(r 1 )r(^)rfe)r(v) F( 3 ) 
rfa+x^+^+v) 



■n::-;-;-:a + n + l,8,|i;P + n + l,A,,^;y + n + l,ri,v; x y z 

-::8 + A,;ri + |^;^ + v: -; -; -; a'b'c 

(3.14) 



Particular Cases: 

Some interesting particular cases of the above results are as follows : 
(i) Taking 8 = oc+l, q = p + 1, £ = y + 1, A, = n = v = nin (3.5), we obtain 

£ £ \' o u c, (l-u) n -'v p (l-v)°-'w 1 '(l-w) n - 1 Y< < "*' ) (xu,yvzw;a,b,c)dudvdw 



r(n)F 



Y, 



(a— n,p— n,y-n) 



(x, y,z;a,b,c) 



(a + l) n (p + l) n (y + l) n n 

which is equivalent to (3.4) 

(ii) Taking 5 = r| = ^ = n+l,X = a, |i=(3, v=yin (3.5), we get 

t t \l u n (l-ur i v I1 (l-v) M w^l-w^ 1 Y i ! a ^ ) (xu,yvzw)dudvdw 

0 _ (n!) 3 



(3.15) 



(«L (PL (y). 



Yr°>(x,y,z) 



(3.16) 



(iii) Replacing 8bya + n+ l— 8, r|byP+n+l — r\, ^byy + n+1— ^, and putting X = 8, \i = r\ and v = ^ in 
(3.5), we get 

J(j *U "U 

_ r( a - 8 + n + i)r(s)r(p - ^ + n + i)t(tiMy - 1 + " + !) r fe) Y^M-r-a (, v , abc ) 

r(a + n + l)r(p + n + l)r(y + n + l) n V ' ' 



Similar particular cases hold for (3.6), (3.7), (3.8), (3.9), (3.10), (3.11) and (3.12). 
(iv) For z = 0, a = b = 2, results (3.4), (3.5), (3.6), (3.7) and (3.9) become 

£ £ x°(r-xr 1 y P (8-yr 1 YM(x,y)dxdy 



(3.17) 



r a + n gP+ n 



Y (a-n, P -n) ^ ^ 



(a + l)„(p + l) n 

£ £ u^fl-u^v^fl-v/- 1 Y^(xu,yv)dudv 

= #ffiffl pl:2;: 
r(8 + ^)r(ri + ia) 



(3.18) 



n: a + n + 1,8; P + n + l,ri; x y 

-: 8 + X ; -q + jj ; 2' 2 

£ £ u^l-urv'll-vr Y n (ra - p) (x(l-u),yv)dudv 

n: a + n + l,X; p + n + l,r|; x y 

-: 8 + X ; r\ + \j. ; 2' 2 

f f ^-'(l-u^V^l-vy 1 - 1 Y^(x(l-u),y(l-v))dudv 

JO Jo 

n: a + n + 1,^; p + n + l,)j,; x y 

-: 8 + X ; r\ + ]j. ; 2' 2 



(3.19) 



_ r(8)r(x)r(ri)r(^) pl:2;2 
r(s+x)r(ri + ia) 



(3.20) 



_ rg^tnW pl:2;2 

r(8 + X)r(ri + ia) 



(3.21) 



IJMER I ISSN: 2249-6645 



www.ijmer.com 



Vol. 4 I Iss. 3 I Mar. 2014 I 4 



A Note on a Three Variables Analogue ofBessel Polynomials 



Results (3.18), (3.19), (3.20) and (3.21) are due to Khan and Ahmad [11]. 
Also, using the integral (see Erdelyi et al. [8], vol. I, p. 14), 

(0+) 



2 i sin n z r(z) = - [ + (- t) z_I e _t dt 



and the fact that 



(1-x-y-z) =2^ 2. L 



r =o s=o j= o r! s! j! 
we can easily derive the following integral representations for Y^™'' 3 ' 1 ^ (x, y ,Z ; a, b, c) : 



(3.22) 



(3.23) 



xu yv zw 
1 + — + — + — 
a b c 



du dv dw 



-r r f (-ur(-vr(-wr n e— 

Jco Jco Jco 

= 8 i (- l) n sin rax sin tcP sin Try r(a + n + l)r(p + n + l)r(y + n + 1) Y n (a ' p ' y) (x, y , z ; a, b, c) 

(3.24) 

(- l) n+1 sin 7ia sin Tip sin Try r(l + a + n)r(l + p + n)r(l + y + n) 



71" 



joo jco jco u _ a _n-l y -p-n-l w -y-n-l g -u-v-w y(o.p,y) 

= (l-x-y-z) n 



ax by cz 

5 5 5 

U V W 



\ 



a, b, c 



du dv dw 

(3.25) 



It is easy to show that 

f (0+) r (0+) f (0+) 



f(u+) p(u+i r(u+) „ 

\ f f u «+n yP+n w Y+n g u+^ H\ 

J— CO J— CO J— 00 



IV. Schlafli's Contour Integral 



1 



xu yv zw 



du dv dw 



J 



v a b c 

= 8 i (- l) n sin rax sin t§ sin Try r(l + a + n)r(l + p + n)r(l + y + n) Y n (a ' p ' y) (x, y , z ; a, b, c) 

(4.1) 

Proof of (4.1) : We have 



(271 



1 f(° +) f( 0+) f ( ° +) u «+n y p + n w y + n e u + v + w yV_ZW_Y 

ri) 3 ^ x ■'- 00 ■'- 00 v a b c y 



du dv dw 



n n— r n— r— s 



v v *v n )r+s+j fx 

"ft ft ft r!s!j! U, 



% s / \ r 



vb. 



1 f(o+) K0+) ko+) 
— u a+n+J v 

j J J-QO J— CO J— CO 



UJ (in 



j p+n+s y+n+r u+v+w 



w T+n + r e u + v + w dudvdw 



(-n) r 



' x w 



+S+J 



n n-r n-r-s 

=y y y 

ft ft ft r!s!j!r(-a-n-j)r(-p-n-s)r(-y-n-r) 
using Hankel's formula (see A. Eerdelyi et al. [8], 1.6 (2)). 



1 1 K0+) 

r(z)~27ciJ = 



e* r z dt 



Finally (4.1) follows from (2.1) after using the result 

r(z)r(l-z) = TT cosecTrz 
for z = 0, a = b = 2, (4.1) reduces to 



(4.2) 
(4.3) 



r(0+) r(0+) ^ a+n v p +n e u+v+w 

J— CO J— GO 



xu yv x 



v 



du dv 



a b 

= - 4 sin rax sin Tip r(l + a + n)r(l + p + n) p) (x, y) 
which is due to Khan and Ahmad [11]. 



(4.4) 
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V. Fractional Integrals 

Let L denote the linear space of (equivalent classes of) complex - valued functions f(x) which are Lebesgue - 
integrable on [0, a], a < oo. For f(x) eL and complex number u. with Rl |i > 0, the Riemann - Liouville 
fractional integral of order u. is defined as (see Prabhakar [13], p. 72) 



r 



f( x ) = — -j-^ f X (x - tf 1 f (t) dt for almost all x e [0, a] 



(5.1) 



(5.2) 



Using the operator F, Prabhakar [14] obtained the following result for Rl u. > 0 and Rl a > -1. 

I" [x« Z: (x; k)]= r r (kn + a + l) ^ Zr (x; k) 
r(kn + a + ia + 1 j 

where Z™ (x; k) is Konhauser's biorthozonal polynomial. 

Khan and Ahmad [11] defined a two variable analogue of (5.1) by means of the following relation : 



and obtained the following result : 

jxj. La+n-x p+„-, Y (a ' p) (x y)l = x ° +ri y P+ " r(a - ^ + n + l)r(p - ia + n + 1) , x 

L * n V y;j r(a + n + l) r(p + n + l) " V 37 

(5.4) 

In an attempt to obtained a result analogous to (5.4) for the polynomial Y ) | a '' 3 ' y ^ (x, y,Z ; a, b, c) we 
first seek a three variable analogue of (5.1). 

A three variable analogue of V 1 may be defined as 

ll " [f ( * y ' z)1 ^ r(x)r(,)r(,) C £ E (x " u) " (y " v) "' (z " w) "' f (u ' v ' w)du dv dw 

(5.5) 

Putting f (x, y,z) - x a+n ~ x y p+n ^ z m Y^ y) (x, y,z;a, b, c) in (5.5), we obtain 

jX^T, = [ x a + n-X yP +n-^ z y + n- n y^.P.y) ^ ^.^ ^ £ )] 

_ x a+ " Y P+n z y+n r(a — A, + n + l)r((3 — ia + n + l)r(y — r| 4- n + 1) {a ^ m) , x 
r(a + n + l)r(p + n + l)r( Y + n + l) n Ix, y,z,a, b, cj 

(5.6) 

VI. Laplace Transform 

In the usual notation the Laplace transform is given by 

L{f(t):s}=J" e st f(t)dt,Rl(s-a)>0 (6.1) 

where f e L (0, R) for every R > 0 and f(t) = 0(e at ), t -> oo. 

Khan and Ahmad [11] introduced a two variable analogue of (6.1) by means of the relation: 



L{f(u,v):s 1 ,s 2 } = £ E [ e s ' u ^f(u,v)dudv 

and established the following results : 



(6.2) 



f- 


2y^ 








V. US 1 




. Sj , s 2 j 



7l 2 S P +n 



sin rax sin7ipr(a + n + l)r(p + n + l) 



(i-x-y) n 



(6.3) 



and 



L 



u a+n y p+n 



+ 



yvs 2 



:s, ,s 2 
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r(g + n + l)r(p + n + l) ( a ,p) , x 

a+n+l B+n+1 n V*» J) 



s l 



(6.4) 



In an attempt to obtain results analogous to (6.3) and (6.4) for Y ] j a '' 3 ' Y '' (x, y,z;a, b, c) we define a 
three variable analogue of (6.1) as follows 

L{f (u, v, w):s 15 s 2 , s 3 } = j" J J e~ s ' U S2V ~ S3W f(u, v, w)dudvdw (6.5) 
Now, we have 



L u 



o-n-l „-p-n~l w -y-"- 1 y(<*.P.Y) I — CZ • 



; a, b, c 



. Sj, s 2 , s 3 



(-l) D+1 7T 3 Sl a+D 4" Sp 



sin rax sin Tip sin Try r(a + n + l)r(p + n + l)r(y + n + 1) 
Similarly, we obtain 



d-x-y) n 



Uu a+n v p+n w y+n 



xus, yvs, zws, 
1 + L + - — - + 



. S|, S 2 ? S3 



r(q + n + l)r^n- t -l)r( T + n + l) ,„,, Tl( 

_ a+n+l _ B+n+1 _y+n+l n V ^ J ' ' ' ' / 



a+n+l B+n+1 y+n+1 
Sj s 2 » 3 



(6.6) 



(6.7) 



VII. Generating Functions 

It is easy to derive the following generating functions for Y i | a ' 3 ' y ^ (x, y,z;a, b, c): 

t ^Y n ^ p - n '^(x,y,z;a,b,c)^e' 
n! 



' xO 



n=0 
00 + n 



< zt vY_1 



V u J 



00 ^.11 

S i T Yi-P--»)(x, y ,z; aj b J c) 
n! 



n=0 



4xt 



1 + Jl- 



4xt 



2yt 
b 



1 + Jl 



4xt 



-p-i 



2zt 
c 



1 + Jl 



4xt 



-Y-l 



(7.1) 



2t 



1 + Jl 



4xt 



(7.2) 



00 ^ n 

S ^-Y i !-P-)(x,y,z;a,b,c) 
^ n! 



I b J 



1 + Jl 



4yt 



2xt 
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Using (3.1), we can also derive the following results : 
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(7.13) 



-U -v -w 



I IJMER I ISSN: 2249-6645 | www.ijmer.com | Vol. 4 | Iss. 3 | Mar. 2014 | 8 | 



A Note on a Three Variables Analogue of Bess el Polynomials 



, xu yv zw 

1 + — + — + 

a b c 



xu yv zw 
1 + — + — + — 
a b c 



. xu yv zw 
1 + — + — + — 
a b c 



00 

X (-l) k ^Y n (-P' 2k - n )(x,y,z;a,b,c) 

k=0 

/•CO (•GO y»00 

= f f | u a+n v 

r(a + n + l)r(p + n + l) Jo Jo Jo 

CO 

X (-l) k ^ 2k Yi- 2k --)(x,y,z;a,b,c) 

k=0 

/•CO |»CC jWGO 

= -? J ua+ " v 

r(a + n + l)r(y + n + l) Jo Jo Jo 

CO 

^ (-l) k l 2k Y n ( 2k -^)(x,y,z;a,b,c) 

k=0 

[ /"CO /•CO /"CO 

~r(p+n+i)r(Y + n + i)J° J° J° u 

CO 

E (-l) k ^ k+1 Y 1 !- p ' 2k+1 -")(x,y,z;a,b,c) 

k=0 

= 1 r r r u -v +n fi+— +^+^ 

r(a + n + l)r(p + n + l) Jo Jo Jo ^ a b c 

CO 

X (-l) k ^ 2k+1 Y;- 2k+1 --)(x,y,z;a,b,c) 

k=0 

| /•CO /•CO /•CO 

~ r(a + n + l)r(y + n + l)J° J° J° " 

CO 

X (-l) k X 2k+1 Y,! 2k+1 - p -)(x,y,z;a,b,c) 

k=0 

| /»CO l"CO 1*00 
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cos A,w du dv dw 
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v 1 
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cos Xv du dv dw 
(7.15) 

cos A,u du dv dw 
(7.16) 

sin A,w du dv dw 
(7.17) 

sin Xv du dv dw 

(7.18) 

sinA,u dudvdw 
(7.19) 



VIII. Double Generating Functions 

The following double generating functions for Y^™'* 3 ' 1 ^ (x, y,z;a, b, c) can easily be derived by using (3.1) 
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IX. Triple Generating Functions 

The following triple generating functions can easily be obtained by using (3.1): 
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X. Bessel Polynomials Of M-Variables 

The Bessel polynomials of m-variables Y^ ai ' a2 ' ' a m) (x^ , X2 , , X m ; aj 



> a 2 ' ' 



be defined as follows: 
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n 
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' a m) 
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(10.1) 



All the results of this paper can be extended for this m-variable Bessel polynomials. The only 
hinderance in their study is the representation of results in hypergeometric functions of m-variables. 
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